The cluster abundance in cosmic string models for structure formation 
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We use the present observed number density of large X-ray clusters to constrain the amplitude of 
matter density perturbations induced by cosmic strings on the scale of 8/i~^Mpc (as), in both open 
cosmologies and flat models with a non-zero cosmological constant. We flnd a slightly lower value of 
(78 than that obtained in the context of primordial Gaussian fluctuations generated during inflation. 
This lower normalization of as results from the mild non-Gaussianity on cluster scales, where the 
one point probability distribution function is well approximated by a x^ distribution. We use our 
estimate of erg to constrain the string linear energy density /i and show that it is consistent with the 
COBE normalization. 



A. Introduction 

Current theories for structure formation can be divided 
into two broad categories: inflation and cosmic defects. 
While inflation predicts primordial and Gaussian fluctu- 
ations (in the simplest inflationary models), topological 
defects [y induce active and non-Gaussian perturbations. 

One of the most important constraints on models of 
structure formation is the the observed abundance of 
galaxy clusters. Although the cluster abundance has 
been widely used to constrain cosmological models with 
primordial Gaussian fluctuations (e.g. [^-||) there have 
been few studies in the context of non-Gaussian pertur- 
bations such as those generated by topological defects. 
This is due to the difficulty of making robust predictions 
in topological defect scenarios, owing to their non-linear 
effects which are difficult to model and require large-scale 
numerical simulations (e.g. Rpl). 

The present work relies on high resolution numerical 
simulations of cosmic string seeded structure formation 
|PHf2|, which are first used to estimate the power spec- 
trum and one point probability distribution functions 
(PDF) of the induced density perturbations (see also 
|jl3Hl6| for different approaches) . We then employ a sim- 
ple generalization of the Press-Schechter formalism p7[ , 
which is suitable for non-Gaussian perturbations with a 
general one point PDF jl^, in order to obtain the ex- 
pected number density of collapsed objects with a mass 
greater that a given threshold. This generalized Press- 
Schechter formalism has been used to constrain the Gaus- 
sianity of the density fluctuations in the Universe ||l9|,Q , 
and has been verified for a particular set of non-Gaussian 
structure formation models including a simplified flat 
space cosmic string model |2^. We finally estimate the 
amplitude of matter density perturbations induced by 
cosmic strings on the standard scale of 8/i~^Mpc, us- 
ing the presently observed number density of large X-ray 
clusters. We do this for cosmic string models in open uni- 
verses without a cosmological constant (50CDM), and 
also in flat universes with a non-zero cosmological con- 
stant (5ACDM) . We use our estimate of erg to constrain 



the string linear energy density /i and show that it is 
consistent with the COBE normalization. 



B. Power spectrum and PDF 

In ref. |1J] we described the results of high-resolution 
numerical simulations of string-induced structure forma- 
tion. The power spectrum of CDM perturbations in- 
duced by long strings can be approximately described 
by 



ScDuiq) = 47rfc37?(fc) oc {O.lqf'-'^l 



where 



p{q) = 3.9- 



2.7 



-0.44 ■ 



(1) 



(2) 



l + (2. 
q = k/T and the shape parameter F is defined as 

r = nl,hexpi^2nl-nynl). (3) 

This fit to our numerical results has an accuracy of bet- 
ter than 10% in the range k = 0.01-100/iMpc~^, provided 
that the baryon to dark matter ratio is relatively small 
(2^2^ + n^/nl^ < 0.3 for 10% accuracy in 7'(fc)). Note, 
however, that there are specified uncertainties in the un- 
derlying numerical results which are significantly larger 
p^ . This result also does not include the contribution 
from small loops which significantly enhance the overall 
power while leaving the overall shape unchanged [ pdj (as 
we will discuss later). A x^ analysis using the obser- 
vational power spectrum reconstructed by Peacock and 
Dodds |£2[ gives a best fit to long string power spectrum 
(0) with F = 0.074 x lO^"-^ at the 95 per cent confidence 
level. For a baryon energy density fls — 0.05, this im- 
plies string models provide a consistent match to obser- 
vations in the acceptable cosmological parameter range 
n'^h = O.Ux 10=^"-^ 

We have used equation (nh to obtain a numerical fit for 
the standard deviation of matter density perturbations at 



the present time as a function of the smoothing scale R 
(in units of h~^Mpc): 



(t{R) = as 



g{8T) ■ 



where 



g(r)=r-'^W, 7(r)=0.5- 



1 + (6/r)' 



(4) 



(5) 



The red-shift dependence of a for arbitrary ilm, ^A is 
described accurately by 



a{R,z) ^a{R) 






(6) 



where 

ff(^m,^A) 





2.5fi„ 




-r!A + (i + an/2)(i + r!A/70) 



(7) 



Here, g{ri^,D,\) gives the suppression of the growth 
of density perturbations relative to that of a critical- 
density- universe [0,p3lp4|, and we can describe the evo- 
lution of the Q-m with red-shift as 



lliYi — \LmyZj 



Kii + ^r 



{i + zni + ni,z-nl) + nl- 



(8) 



Using the simulation results for string-induced matter 
perturbations, we also found a reasonable fit to the pos- 
itive side of one point PDF as follows: 



X^(?^) — n 



2n 



(9) 



where x^ ("■) is the standard x^ variant with the number 
of degrees of freedom n, given as 



n^0A{RT + 2) 



3.6 



(10) 



Here R is again the smoothing radius in units of h ^Mpc 
and the fit is accurate within 10% error for all scales. 



C. Modified Press Schechter formalism 

The Press-Schechter formalism [O relates the mass 
fraction of collapsed objects whose masses are larger 
than some given threshold M, with the fraction of space 
in which the evolved linear density field exceeds some 
threshold 6c- It has been extensively tested with great 
success against N-body simulations in the context of pri- 
mordial Gaussian fluctuations, allowing for the compu- 
tation of the number density of clusters within a given 
background cosmology. Here we used an extension of 



the Press-Schechter formalism for non-Gaussian pertur- 
bations proposed by Chiu, Ostriker and Strauss |1^ and 
verified for a particular set of non-Gaussian structure for- 
mation models in ref . pl| . The fraction of the total mass 
within collapsed objects larger than a given mass M is 
given by: 



an(>M(i?,z),z) 

rimiz) 



Ir'x PR{v)dv, 



(11) 



<5. 



where M is the cluster mass defined by M = 47r_R'^/9b/3, 
fn, = l//o PR{v)dv, Pr,{v) is the one point PDF for a 
given smoothing radius R, v is the number of standard 
deviations from mean density, (5* = 6c/ ^{R, z) and 6c = 
1.7 ± 0.2 (95% confidence interval) assuming spherical 
collapse p5| . 

To obtain the number density of clusters in a mass 
interval dM about M at a redshift z, we differentiate the 
Press-Schechter formula (O) to obtain (c.f. ref. M) 



n{M, z)dM 
~fR^ 



da{R,z) 



M <7^(R,z) dM 



Pn 



<t(R./. 



dM. 



(12) 



In the derivation of equation (03) we have ignored other 
terms resulting from the R dependence of P/j and Je.- 
This is a mathematically motivated approximation in the 
regime of mild non-Gaussianity, and we have verified in 
our case that it gives rise to at most an extra 2% error in 
the final estimate of erg . Substituting (H) into (|l^ ) gives 



n(M,z)dM = /^3^^ 



f o.76(flr)°-^iog(it'r) 
\ [{RTy>--i +i.7i\-^ 



2a^(R,z) 

+ -i{RT)\dM. 



(13) 



Lacey and Cole constructed a merging history for dark 
matter halos based on the excursion set approach and 
obtained an analytical expression for the probability that 
a galaxy cluster with present virial mass M would have 
formed at some redshift z ||2^. The probability that a 
galaxy cluster with present virial mass M would have 
formed at a given red-shift z is given by: 



p{z) ^p{w{z)) 



Aw{z) 



where 



p{w{z)) 



w{z) 



2w{z){F-'-l)eTic(^ 

y|(F-i-2)exp(-^ 
_ 4(g(M,0)/g(M,z)-l) 



(14) 

(15) 
(16) 



and F — 0.75 ± 0.15 (95% confidence interval) is the 
fraction of the cluster mass assembled by a red-shift z 
pTj . Wc note that this result was derived in the context 
of primordial Gaussian fluctuations and must ultimately 



be verified using N-body simulations. Although it may 
not be valid for generic non-Gaussian models, we still 
expect this to be a good approximation in the context 
of the cosmic string scenario for structure formation as 
the departures from Gaussianity on clusters scales are 
relatively small |l^ . 



D. The Mass-temperature relation 

In order to use the generalized Press-Schechter formal- 
ism to determine the abundance of X-ray clusters with a 
given temperature, we need to relate the X-ray tempera- 
ture of a cluster with its virial mass. Here, we use the re- 
sults of ref . y for the normalized virial mass-temperature 
relation (modified from H|): 



Mv 



(1.23 ±0.33) X IQis 

2{2-r,)(4-ny 



m 



kaT 



1/2 



1.67 

1 + Zt 



64-56?)+24r;2-7j;3+77<i 7.5koV 



3/2 



(17) 



h-^M. 
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where zt is the turnaround redshift, r2,nt = ^m(^t), and 
32 



V = Vizt) 



"A"mt 



b{n) 



9^2f70^(l+^^)3' 

0.76-0.25f7 (OCDM), 
0.73-0.231} (ACDM). 



(18) 
(19) 



For a given red-shift of cluster collapse Zc, the turnaround 
redshift zt is easily obtained using the fact that 2t(zt) = 

i(zc). 

Hence we can now estimate the present comoving num- 
ber density of galaxy clusters per temperature interval 
d(fcBr) with a mean X-ray temperature k^T which were 
formed at a given redshift z as; 



tiTik^T, z)A{k^T)dz 



3M 
2kBT 



n{M, 0)p{z)d{kBT)dz. 



(20) 



The present abundance of X-ray clusters with a temper- 
ature fceT greater than 6.2keV can be estimated by in- 
tegrating equation (^0|) from z = Zc = to infinity. A 
comparison between the observed cluster abundance and 
its theoretical prediction will give an estimate of ag. 

We will use the observation for the number density 
of galaxy clusters at z = 0.05 with X-ray temperature 
exceeding 6.2keV, given by Viana and Liddle B, based 
on the dataset presented by Henry and Arnaud ]28[, and 
updated by Henry [g9| : 

N{> 6.2keV, 0.05) = 1.53 x 10"^"-^^/i^Mpc"^ (21) 

The uncertainty in ( pl| ) is the 1-sigma interval, and these 
results have taken into account the effect of temperature 
measurement errors. The reasons for concentrating on 
galaxy clusters with temperature larger than 6.2keV has 
been extensively discussed by Viana and Liddle |4J . 




FIG. 1. The ug in cosmic string scenarios normalized to 
cluster abundance. The upper panel is the iSOCDM models, 
while the lower panel is the iSACDM models. 



E. Results and discussion 



By comparing IGW) with the result integrated from 
(PQ), we obtain the observationally constrained ug as 
plotted in figure m The overall error in the value of 
as was estimated by Monte Carlo simulations over 10^ 
realizations, treating the intrinsic uncertainties in F, 
iV(> 6.2keV, 0.05) as lognormal, and those in 5c, My, 
F as Gaussian. An accurate numerical fit to this result 



IS 



0-8 



0.44 rj' 
0.44 rj' 



-0.45+0. isn" 



-0.6+0.3O„ 



(50CDM), 

(5ACDM). 



(22) 



The 95% confidence limits are ±32% in the 50CDM case, 
and +35% and -32% in the 5ACDM case. We note that 
the overall shape of erg in figure Olis more or less the same 
as that for inflationary models^], while the amplitude 
here is about 10-20% lower than standard inflation. The 
emergence of the same shape from these very different 
scenarios is due to the fact that both the string-induced 
and inflationary power spectra used in the calculation of 
cluster abundance are constrained by the same observa- 
tion p2| and have roughly the same shape within the 
scales of interest (though with quite different choices of 
F). We have also verifled that the lower normalization of 
(78 is due to a slightly larger right-hand side tail of the 




FIG. 2. Bg/j. as a function of fi^ and rj^v- The hatched 
(50CDM) and shaded (5ACDM) areas are within a 95% con- 
fidence level. 

PDF in the cosmic string case, i.e. a high-density region 
then becomes consistent with a smaller as- 

The value of as in the context of the cosmic string 
model for structure formation was also investigated by 
Bruck [Sy]. In his work he assumed the one point PDF 
to be independent of scale by taking the distribution at 
the non-Gaussian scale (~ 1.5(r2^/i^)"^Mpc 0) to be 
valid up to scales relevant for the cluster abundance cal- 
culation. Although this assumption can give the right 
qualitative results, for small values of fij^/i, the R de- 
pendence of the one point PDF needs to be taken into 
account in order to obtain more accurate results. This 
improvement has been incorporated in our work, which 
also took into consideration the merging history of dark 
matter halos. 

To see how cluster abundances constrain the string en- 
ergy density per unit length /z, we can define a string 
"bias" parameter, 



BgM 



" ^1) 



as{^\ 



.Ml) 



.(S) 



i^'n.^1) 



(23) 



.(•5) 



where trg is the cluster result (G3) and a^ is directly in- 
ferred from the long string-induced power spectrum p), 
given the assumptions that fig = 0.05, ri°-,ft, = 0.14, and 
G/i6 = G/i X 10^ = 1 for nl-^ = 1 with A = (with 
the open and A-parameter dependence for GfXe given by 
Gfi6 oc flH, °-^ for 50CDM and G/ie ex n°^~°-°^ for 
5ACDM Q). 

We can observe in figure for the flat A-model that 
BG^{n'^,nl) ranges from 2.2±0.8 to 3.9±1.3 (95% con- 
fidence level) for different choices of 17j^ and 17^. These 
values of Sg^ are actually the cluster normalized values 
for G^le at 1^;^ = 1 with A = 0. For small fij^ w 0.2-0.3, 
the result is slightly higher than the previously COBE- 
normalized G^e = 1.7 [ pl[ , but it is certainly consis- 
tent within the uncertainties. For large fi^, the result 



appears to be inconsistent with the COBE constraint. 
However, as recently discussed in ||ll|, the inclusion of 
perturbations induced by cosmic string loops can boost 
the string-induced power spectrum by a factor of about 
two under reasonable assumptions. This may serve to re- 
move even this apparent discrepancy in G/iq between the 
COBE and cluster abundance constraints. We also note 
that the exact contribution of the background of gravita- 
tional radiation emitted by cosmic string loops remains 
a significant uncertainty 132] . 



F. Conclusion 

We have constrained the amplitude of matter den- 
sity perturbations induced by cosmic strings on the scale 
of 8/i^^Mpc in both open cosmologies and flat models 
with a non-zero cosmological constant, using the cur- 
rently observed number density of large X-ray clusters. 
Because string-seeded matter perturbations are mildly 
non-Gaussian on cluster scales, we obtained a slightly 
lower normalization of as than that found for cosmolog- 
ical models with primordial Gaussian fluctuations. We 
used the calculated as to constrain the string linear en- 
ergy density /i, which we found to be consistent with the 
latest COBE normalization when current uncertainties 
in the normalization are taken into account. 
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